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ABSTRACT 

We investigate thermal instabihty in an electron-ion magnetized plasma relevant to 
galaxy clusters, solar corona, and other two-component astrophysical objects. We apply 
the multicomponent plasma approach when the dynamics of all the species are considered 
^ ■ separately through electric field perturbations. General expressions for perturbations 
Q ■ obtained in this paper can be applied for a wide range of multicomponent astrophysical 
^ ■ and laboratory plasmas also containing the neutrals, dust grains, and other species. We 
assume that background temperatures of electrons and ions are different and include the 
-T— , . energy exchange in thermal equations. We take into account the dependence of collision 
^ ! frequency on density and temperature perturbations. The cooling-heating functions are 

taken as different ones for electrons and ions. As a specific case, we consider a condensation 
mode of thermal instability of long-wavelength perturbations when the dynamical time is 
smaller than a time during which the particles cover the wavelength along the magnetic 
field due to thermal velocity. We derive a general dispersion relation taking into account the 
effects mentioned above and obtain simple expressions for growth rates in limiting cases. 
Perturbations are shown to have an electromagnetic nature. We find that at conditions 
under consideration transverse scale sizes of unstable perturbations can have a wide 
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spectrum relatively to longitudinal scale sizes and, in particular, form very thin filaments. 
The results obtained can be useful for interpretation of observations of dense cold regions 
in astrophysical objects. 

Key words: conduction - galaxies: clusters: general - instabihties - magnetic fields - 
plasmas -waves 
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1. INTRODUCTION 

The thermal instabihty leads to formation of regions with larger densities and lower 
temperatures than that in the surrounding medium (Parker 1953; Field 1965). Beginning 
from the classical paper by Field (1965), this instability was studied for both astrophysical 
objects (for reviews see, e.g., Vazquez- Semadeni et al. 2003; Elmegreen & Scalo 2004; 
Cox 2005; Heiles & Crutcher 2005) and plasma physics applications (e.g., Meerson 1996). 
Majority of papers were devoted to thermal instability in the interstellar medium (ISM; 
e.g.. Field 1965; Burkert & Lin 2000; Hennebelle & Perault 2000; Koyama & Inutsuka 
2002; Kritsuk & Norman 2002; Sanchez-Salcedo et al. 2002; Audit & Hennebelle 2005; 
Stiele et al. 2006; Vazquez-Semadeni et al. 2006; Fukue & Kamaya 2007; Inoue & Inutsuka 
2008; Shadmehri et al. 2010). Solar prominences are supposed to be formed as a result of 
thermal instability (e.g.. Field 1965; Nakagawa, 1970; Heyvaerts 1974; Mason & Bessey 
1983; Karpen et al. 1989). In galaxy clusters, this instability, including the presence of the 
magnetic field, was studied in (e.g.. Field 1965; Loewenstein 1990; Balbus 1991; Bogdanovic 
et al. 2009; Parrish et al. 2009; Sharma et al. 2010). The nonhnear stage of thermal 
instability resulting in formation of nonlinear cool structures was investigated in the ISM 
(e.g., Trevisan & Ibaiiez 2000; Sanchez-Salcedo et al. 2002; Yatou & Toh 2009) and solar 
corona (Mason & Bessey 1983; Karpen et al. 1989; Trevisan & Ibaiiez 2000). 

In papers studying thermal instability in astrophysical objects with the magnetic field, 
the one-fluid ideal MHD is generally used. The two-fluid model of the ideal MHD has been 
treated, e.g., by Fukue & Kamaya (2007) and Inoue & Inutsuka (2008). The non-ideal 
effects in the magnetic induction equation have been considered by several authors (e.g., 
Heyvaerts 1974; Stiele et al. 2006; Shadmehri et al. 2010). 

For astrophysical media consisting of many kinds of species (electrons, ions, dust 
grains, neutrals, and so on), the multicomponent approach considering the dynamics of 
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each species separately is an adequate method of investigation (e.g., Nekrasov 2009a, 
2009b, 2009c). The thermal instability in multicomponent media has been studied by 
Kopp et al. (1997), Pandey & Krishan (2001), Pandey et al. (2003), Shukla & Sandberg 
(2003), Kopp & Shchekinov (2007). Analytical investigation of thermal instabihty in 
multicomponent magnetized media with such physical effects as collisions between different 
species, ionization and recombination, dust charge dynamics, gravity, self-gravity, and so 
on is a sufficiently difficult problem. Therefore, one usually treats simplified models such 
as, for example, potential perturbations in nonmagnetized (Kopp et al. 1997; Pandey & 
Krishan 2001; Ibanez & Shchekinov 2002; Pandey et al. 2003; Shukla & Sandberg 2003; 
Kopp & Shchekinov 2007) and magnetized (Kopp et al. 1997; Shukla & Sandberg 2003) 
plasmas. 

When studying thermal instability, one usually does not take into account an energy 
exchange between species in thermal equations. It may be done at a weak or strong 
coUisional coupling of species. However, an intermediate case can in general also occur. 
The inclusion of this effect results in considerable analytical complications (e.g., Birk 2000; 
Birk & Wiechen 2001). The absence of thermodynamical equilibrium is an additional factor 
complicating a problem. However, different temperatures of species can be observed, for 
example, in galaxy clusters (Markevitch et al. 1996; Fox & Loeb 1997; Ettori & Fabian 
1998; Takizawa 1998). Therefore, this effect needs also to be taken into consideration. 
When background temperatures of species are different, it is necessary to take into account 
the perturbation of the energy exchange frequency which depends on the number density 
and temperature. 

In this paper, thermal instability in the electron-ion magnetized plasma relevant to 

galaxy clusters, solar corona, and other two-component astrophysical objects is investigated. 
We apply the multicomponent plasma approach when the dynamics of all the species 
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are considered separately through electric field perturbations (the E-approach; see, e.g., 
Nekrasov 2009a, 2009b, 2009c; Nekrasov & Shadmehri 2010, 2011). General expressions 
obtained in this paper can be applied for a wide range of astrophysical and laboratory 
plasmas also containing the neutrals and dust grains. We assume that background 
temperatures of electrons and ions are different and include the energy exchange in thermal 
equations. We take into account the dependence of energy exchange collision frequency on 
density and temperature perturbations. The cooling-heating functions are also considered 
for both electrons and ions. We do not include ionization and recombination effects and 
the gravity. Expressions for electron and ion perturbations are obtained in the general form 
which can be used for other species. As a specific case, we here treat a condensation mode 
of thermal instability of perturbations elongated enough along the background magnetic 
field. In this case, the dynamical time is smaller than a time during of which the particles 
cover the longitudinal wavelength due to their thermal velocity. The opposite, fast sound 
speed limit, is considered in (Nekrasov 2011). We derive the general dispersion relation, 
taking into account the effects mentioned above, and discuss the limiting cases. 

The paper is organized in the following manner. In Section 2, we give fundamental 
equations used in this paper. An equilibrium state is considered in Section 3. General 
equations for temperature perturbations are obtained in Section 4. In Section 5, 
equations for components of velocity perturbations are given in the fast dynamical regime. 
Components of perturbed current are calculated in Section 6. These components for the 
simplified collision contribution are given in Section 7. In Section 8, we derive the dispersion 
relation. Its limiting cases are considered in Section 9. We discuss the obtained results in 
Section 10. The possible astrophysical implications are considered in Section 11. Summing 
up of main points is given in Section 12. 
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2. BASIC EQUATIONS 



The fundamental equations we use are 



^'■+v,.Vv, = -^ + F,+:^v,xB, (1) 



the equation of motion, 

^ + V • n,v, = 0, (2) 

the continuity equation, 

dT- 1 

--^ + • VTi + (7 - 1) r,V • = - (7 - 1) -d {m, Ti) + i/^ (ne, Te) (Te - T^) (3) 

OT Hi 

and 

^+Ve-VTe+(7 - 1) TeV-Ve = - (7 - 1) - V ■ q,-(7 - 1) -^e (ne, Te)-z/|, (n„ T^) (T^ - T,; 

(4) 

are temperature equations for ions and electrons. In Equations (1) and (2), the index 
j —i,e denotes ions and electrons, respectively. The value Fj in Equation (1) is given by 



F. = - ly,^ (v, - Ve) , (5) 

rrii 

Fe = — E - Uei (Ve - Vj) . 

me 

Other notations in Equations (l)-(5) are the following: qj and rrij are the charge and mass 
of species j = i, e, \j is the hydrodynamic velocity, rij is the number density, Pj = UjTj is 
the thermal pressure, Tj is the temperature in the energy units, i/ie {uei) is the collision 
frequency of ions (electrons) with electrons (ions), uf^{ne,Te) — 2uie {vl^{ni,Te)) is the 
frequency of the thermal energy exchange between ions (electrons) and electrons (ions) 
(Braginskii 1965), nivf^ {ne,Tg) — Ue^li {ni,Te), 7 is the ratio of specific heats, E and B 
are the electric and magnetic fields, and c is the speed of light in vacuum. The value qe 



-7- 



in Equation (4) is the electron heat flux associated with the thermal motion in the system 
of coordinates where the electron gas is at rest as a whole (Braginskii 1965). As for the 
latter, we will consider a weakly coUisional plasma when the electron Larmor radius is much 
smaller than the electron coUisional mean free path. In this case, the electron heat flux is 
mainly directed along the magnetic field, 

qe = -Xeb(b-V)re, (6) 

where Xe is the electron thermal conductivity coefficient and b = B/i? is the unit vector 
along the magnetic field. In other respects, a relation between cyclotron and collision 
frequencies of species stays arbitrary in general expressions considered below. We only take 
into account the electron heat flux (6) because the corresponding ion thermal conductivity 
is considerably smaller (Braginskii 1965). We also assume that the heat flux in equilibrium 
is absent. The cooling and heating of plasma species in Equations (3) and (4) are 
described by function Cj{nj,Tj) — rijAj (Tj) — n^Tj, where Aj and Tj are the coohng 
and heating functions, respectively. The form of this function differs from the usually 
used cooling-heating function £, beginning from the classic paper by Field (1965). Both 
functions are connected with each other via equality Cj {nj,Tj) — mjnj£j. Our choice is 
analogous to that as in Begelman & Zweibel (1994), Pandey & Krishan (2001), Shukla & 
Sandberg (2003), Bogdanovic et al. (2009), Parrish et al. (2009). The function A^- (T^) can 
be found, for example, in Tozzi & Norman (2001). 

Electromagnetic equations are Faraday's 



and Ampere's 



VXE = --- (7) 



V X B =— j (8) 

c 



laws, where j = Ylij We consider wave processes with typical timescales much larger 

than the time the light spends to cover the wavelength of perturbations. In this case, one 
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can neglect the displacement current in Equation (8) that results in quasi-neutrality for 
both electromagnetic and purely electrostatic perturbations. The magnetic field B includes 
the background magnetic field Bq, the magnetic field Bocur of the background electric 
current (when is present), and the perturbed magnetic field. 

For generality, we assume in the meanwhile that 7^ ng, having in mind that some 
expressions obtained below can be applied for multicomponent plasmas. 

3. EQUILIBRIUM STATE 

At first, we will consider an equilibrium state. We assume that the background 
fiow (average) velocities of species are absent. We do not here involve an equilibrium 
inhomogeneity. Then, thermal equations (3) and (4) in equilibrium take the form 

(7 - 1) —A {riio, T,o) - lyU^eO, T,o) {T,o - T,o) = 0, (9) 
_ 1) J-C, (rieo, Teo) + v% (n,o, T,^) (Teo - Tio) = 0, 

^eO 

where the subscript denotes equilibrium values. 

4. LINEAR EQUATIONS FOR TEMPERATURE PERTURBATIONS 

We now consider Equations (3) and (4) in the linear approximation. Applying the 
operator d/dt to Equation (3) and using for ions Equation (2) to exclude the number 
density perturbation and Equation (9), we find 



(10) 
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where and below the subscript 1 denotes perturbed values. The operators and notations 
introduced in Equation (10) are as follows: 



9 _ ^ \ d 



d 

D2i — {^Tie + ^ie) 

d 



Cii — Tjo 

C2i = ^lie {TgQ — Tjo) . 



(7 - 1) ^ + ^ni 



{TeO — Tio) 



iO 



(11) 



Analogously, we obtain for electrons 



DleTei - -026^11 = C'leV ■ Vgi + C2eV • Vji, 



(12) 



where 



d \ d 

Die — ( TTT + + ^Te + ^Tei + ^ei ) 



— ^ei 



dt ■ 

d_ 



Cie — 
C2e = ^ei {TgQ — Tjo) . 



(7 - Ij ^ + i W H i^e 



(13) 



In notations (11) and (13), we have introduced the following frequencies: 

XeO d'^_ 
2 ' 

dCi {nio,Tio) 



fi^ = -(7-l) 

^Te = (7 - 1 



dCe (^ecTeo) 



neo<9T, 



,^^ri = (7-l) 



eO 



niodTio 



^ne = (7 - 1) TfT-^ , = (7 " 1) 

^ei = ^'ei (^iO) ^eo) , ^^ie = ^ie (^eO, T'eo) , 



Tioduio 



^Tei — 



(.-t eO — -t iOj , "Tie — ^77; (.-ieO " -tiOj 



eO 



eO 



(14) 



- 10 - 



We assume that the background magnetic field Bq is directed along the z-axis. In notations 
(11) and (13), we have used the equilibrium state and the number density dependence 
of u^- {nio,Teo) ~ TiiQ and uf^ (^eo^^^eo) ~ ^eo- We see from Equations (10) and (12) that 
temperature perturbations are connected with a velocity divergence. Solutions for Tgi and 
Til are given by 

DTei = Gi V-Vei + GaV-va, (15) 
L>ra = G'3V-Vei + G'4V-Va, (16) 
where the following notations are introduced: 

D = {DuDi, - D2iD2e) , (17) 

Gi = {D\iC\e — D2eC2i) , 

G2 — {DiiC2e + D2eCii) , 

G3 = {D2iCie — DieC2i) , 

G4 — {DieCii + D2iC2e) ■ 

To find the temperature perturbation T!,i, we must have expressions for V-Vji. General 
equations for the velocity v^i and V • v^i are derived in the Appendix, where expressions 
for D and G/, / = 1,2,3,4, are also given. In their general form, the components of Vji 
are very complex. Therefore to proceed further analytically, we here restrict ourselves to a 
limiting case in which the dynamical time {d/dt)~^ is short in comparison with a time the 
thermal particles need to cover the wavelength along the magnetic field. Some additional 
simplifying conditions which are satisfied in magnetized plasmas are also used. 

We note that the opposite case when the dynamical frequency is smaller then the 
corresponding sound frequency has been considered in (Nekrasov 2011). The general 
expressions are the same in the last and this papers. However for convenience of reading, 
we keep their here. 



- 11 - 



5. SPECIFIC CASE: ft > ^ (4e + 4J 

Equations (A27), (A28) and (A30) are written in their general form, which allows us 
to consider different simplified specific cases corresponding to real astrophysical conditions. 
We further proceed with sufficiently fast perturbations such that 

This condition is opposite to the one for the fast sound regime (Nekrasov 2011) and 
corresponds to the long- wavelength perturbations along the magnetic field. At the same 
time, we assume that ^ d^/dt^ for magnetized plasma. Other condition is a common 
one for hydromagnetic description, i.e. 

when the Larmor radius of species is much smaller than the transverse wavelength of 
perturbations. The square of the velocity v"^^ has the following estimation: 

^2 ^ ^eO {§i + ^ie) + TiQQ.^i ^^^^ 



'^i {it + + ^ei) 

Expressions (20) arc given in the approximate form to unite two cases, ^ > (or <)f2je,ej- 
We note that operators {d/dt)~^ and {d/dr)~^ in expressions (18)-(20) and corresponding 
expressions below denote typical dynamical times and wavelengths of perturbations. Under 
conditions (18) and (19) and using notations (AlO), (A30), and (A31), we find equations 
for Pi,ei(see Equations (A27) and (A28)): 

/I d'^Fiix dFiiA f 1 a^Fei^ , , 

^a = Aj-— —^ + —^ —^ + —^ , (21) 



Uci dydt dz J \ Uce dydt dz 

p \ [ ^ ^^-^elx 9Feiz\ f 1 d'^Filx 9Fiiz 1 (22) 

^ ' ojce dydt dz J \ Ud dydt dz ' ' 
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where the following notations are introduced: 









G3 


nii 




DlTLi 












me 


[dtj 


1 i^me'^^ 


Drrie 



(23) 



Equations (21) and (22) have a symmetric form relatively to changing the index i by e and 
vice versa. Estimations of \j and jij are followed from expressions (A16), (A18), (A19), and 
(A29). 



5.1. EQUATIONS FOR COMPONENTS OF VELOCITIES v^ei 
We now obtain equations for components of velocities Vj^ei) using Equations (21) and 



(22). 



5.1.1. Equations for Vi^eiy 

Prom Equations (A3), (21), and (22), we find, using notations (A6), 

_ _ 1 P 1 d^Fiix 1 dFay 1 d^Fay 

^ ilx ^ elx 



ul^dy'^dt \Uei iOc 



1 92 

ul dydz 



1_ 1 d'^Fei^ J_dFeiy _ J_d^Feiy 



ce ^ce ^ce ^ce 

3 . . 



1 d-" f _^ 



ce 

■^2 



0(1 c;±u- 



1 8^ 

r Q Q ('^e-^elz — fJ'eFilz) ■ 
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The terms proportional to are needed in equations for Vi^eix for obtaining terms ~ '^cj'- 
We see that these solutions are obtained one from another by changing i •<-> e. 

5.1.2. Equations for Vi^eix 
Equations for Ui,eia; are easily found from Equation (A2) by using Equations (24) and 



(25): 



Vix. = —Fay + — ^ - —^7? 26 

^ ilx ^ elx 



ce 



1 d"^ fdY^ 



r\ o \ tlx 

I fd_ 

LOci dydz \ dt 

II dFeix 1 d^Fe 

Velx — ^ely H 5 5 H72~ ^^Tj 

ujce ujI dt ul^ dt^ 

1 d'' [ K 

^ e\x ^ ilx 



1 d"^ f d 



1 /_a 

Wee Sy^z \ ^ 

5. i . 5. Equations for Vi^^iz 

From Equations (A7), (21), and (22), we obtain, using notations (A6), 

9\iz ^ dFiu _ / ^Fax - —Feia) (28) 

dt'^ dt dydzdt \ojci ^ ^ uJce ^ ^ ) 

aH.._3F,,. /A^^^_^_^ X ^^^^ 



dt'^ dt dydzdt \ujce 

+ {KFelz — ^eFilz) 
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6. COMPONENTS OF CURRENT 



We now find components of the finear current ji = ^^?j%oVji. It is convenient to 

j 

consider the value At: [d j dt)~^ In this case, we obtain dimensionless coefficients by Ei. 
We further consider the electron- ion plasma in which rieo = fiio — no, Qe — —qi- In our 
calculations, we will use an equality me^ei = fniVi^ and the relation mjFji = — meFgi. Prom 
Equations (24)-(29) and using notations (5) in the linear approximation, we find 



(30) 



bxx {Vilx - Velx) + bxy {Vily - Vely) - bxz {Vilz " Velz) , 



( — ) jly = ttyxElx + dyyEly + ttyzElz 



(31) 



byx {'^ilx ^elx) byy {Vily ^ely) byz (fjlz ^elz) ; 



47r 



jlz — dzxEix + ClzzEiz 



(32) 



- bzx {Vilx - Velx) - bzz {Vilz " Veu) , 



Here, the following notations are introduced: 



= ^ <( 1 [(Ae - He) rrie + (Aj - //j) rrii 



dW d 



col 



rrii 



dy"^ \dt 



~ col dt 



pi 



XiTTli - lleTTle XeTUe - fiiTTli 



rrii 



nip 



( d 



-2 



^yx — 



^pi \ \ d 1 



nii 



(Ai - \Xi) mi (Ae - lie) m, 



COr. 



CO, 



dyz — 



pi 



1 / X rrii 
— Aj + jii — 



+ 



Ae + A*, 



OJr. 



nii 



dydz \dt 



f d 



dydz \dt 



[(Ae - He) - [K - Hi)\ ^ _ n „ 

UJri 



(33) 



- 15 - 



where Upi — {ATrriiQqf /rrii)^^^ is the ion plasma frequency. 



7. SIMPLIFICATION OF COLLISION CONTRIBUTION 

Relationship between uuce and or oud and (that is the same) can be arbitrary 
in Equations (30)-(32) (except of that in the thermal conduction). We further proceed 
by taking into account that d/dt ^ Ud- In this case, we can neglect collisional terms 
proportional to b^y and by^ (see notations (33)). However, a system of Equations (30)- (32) 
stays sufficiently complex to find ji through Ei. Therefore, we further consider the case in 
which the following condition is satisfied: 

1>^TT7<1 K (Ae - /ie) + mi (A, - /X,)] — — ) (34) 



at I rrii oy^ \ot 

or on the order of magnitude 

f d 



1»^^ 



where = (Teo + T^o) /mi. It is obvious that this inequality can easily be realized in 
magnetized plasma. Under condition (34), we can neglect the terms ~ bxx and bzx in 
Equations (30) and (32). In the case » Uigd/dt, the term ~ byy in Equation (31) can 
also be omitted. Thus, a system of Equations (30)-(32) takes the form. 



47r ( — j jix — SxxElx — £xyEly + ExzEiz, (35) 



ATT ( — J jly = EyxElx + £yyEly + EyzElz, (36) 

47r j jiz = SzxEix + SzzEiz- (37) 
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The following notations are here introduced: 



where 



fip OjyzOjzx ^yz 



ujidtii + dY^"' ^'"'^"^ (i + d)' 



pi 



"pi 

Parameter d defines the collisionless, o? <C 1, and collisional, o? 3> 1, regimes. Below, we 
derive the dispersion relation. 



8. DISPERSION RELATION AND ELECTRIC FIELD POLARIZATION 

We further consider Equations (35)-(37) in the Fourier-representation, assuming that 
perturbations have the form exp (ik • r—iut). Then using Equations (7) and (8), we obtain 
the following system of equations: 

(n^ — Exx) Eixk + ^xyEiyk — SxzEizk — 0, (40) 

—EyxEixk + {nl. ~ ^yy) Eiyk — {riyHz + Syz) Ei^k = 0, 
—SzxEixk — nyUzEiyk + (riy — Ezz) Eizk — 0, 

where Ei^ is the Fourier-image of the electric field perturbation, n = kc/u;. The index k by 
Eifc is equal to = {k,a;}. For the Fourier-images of operators £ij and c/, we keep the same 
notations. In general, wc sec that the longitudinal electric field Ei^k ~ Eix^yk inevitably 
arises when ky ^ and riy — ezz ^ 0- The dispersion relation can be found by setting the 
determinant of the system (40) equal to zero. 
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We will consider the case in which 

Szz > nl. (41) 

In the coUisionless regime (d <^ 1), this inequality denotes that the transverse wavelength of 
perturbations is much larger than the electron skin-depth. We assume that condition (41) 
is also satisfied in the coUisional regime (c? ^ 1). Expressing the electric field through 
Eix,yk from the third Equation (40) and substituting it into two other equations, we find, 
using expressions (33) and (38), that contribution of the longitudinal electric field Ei^k in 
the case (18) is neghgible. For an estimation of values (33), we take 

(A,-/x,)m,~(r,o + r,o)(^) (42) 

(see expressions (23)). Prom Equations (33) and (38), it is easy to see that Syxj^xx ^ oj/ud 
and Exy/Syy ~ oj/uJci- Since uj"^ <^ cu^j, the contribution of term ExySyx into dispersion 
relation is small as compared with that of term Sxx^yy Then using an estimation (42), we 
obtain the simple dispersion relation 

{nl - Syy) (n^ - Sxx) ^ 0. (43) 

The first factor of Equation (43) describes the Alfven wave cu"^ = klc\, where 
ca — Bq/ (Anmiriio)^^^ is the Alfven velocity. The polarization of this wave is 
El = (0, Ely, 0). The second factor of Equation (43), 

- Sxx = 0, (44) 

describes the magnetosonic kind of perturbations with polarization Ei = (£'ij;,0, 0). This 
perturbation is purely the electromagnetic one. 



9. SOLUTION OF DISPERSION RELATION (44) 
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The parameter d defined by Equation (39) is equal to d — iuei/uj. li d <^ 1, then 
the coUisional term in e^^ is imaginary one and of the order of {u^i/uj) {k^v^^/u'^) <^ 1 in 
comparison with the term a^x- In the case d ^ 1, this collisional term becomes the real one 
and is ~ {klv^^/uj^) <^ 1 in comparison with axx- Therefore, the collisional term in Exx is 
not taken into consideration in Equation (44). Then, the dispersion relation takes the form, 

00^ = ec\ + [D {Teo + Tio) + iuj (d + + + G4)] . (45) 

The values D and Cj, i — 1,2, 3, 4, are given by expressions (A16)-(A20), where d/dt should 
be replaced by —iou. Calculating expression in the square brackets in Equation (45), we 
obtain 

u^'-k^cl + kl^. (46) 

Here, 

R = TeO (-i-fU + n^ + Qre - ^ne) {-i^ + ^Ti + '^^ie) (47) 
+ Tjo {-i'^U + VtTi - flni) {-i^ + + ^Te + 2f^ei) 
+ (TeO - Tio) (-ioj + n^ + Qtc) ^ie + (^iO - Teo) {-lU + ^Ti) ^ei 

+ TeO {-i^ + ^Ti) ^Tei + TeO [-i (7 - 1) - fine] ^Tie + ^iO (-^7'^ + ^Ti - ^ni) ^Tei 

and 

V = [{-iu + + ^Te) {-ioO + ^Ti + fiie) + {-ioO + ^Ti) {^ei + ^Tei)] , (48) 

where = (7 — 1) (Xeo/^eo) k'^. All the values Q are defined by a system (14). We see that 
the first four terms on the right hand-side of expression (47) are symmetric ones relatively to 
contribution of electrons and ions. The last three terms are connected with perturbation of 
collision frequency {rii^eo, ^e) because of the electron temperature perturbation. All the 
terms proportional to fiei,ie and ^lTei,ie in expressions (47) and (48) are connected with the 
energy exchange in thermal equations (3) and (4). We see that this effect in a general case 
when 7^ T'io results in considerable modification (complication) of dispersion relation. 
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Nevertheless, it must be taken into account because the absence of thermal equilibrium 
between electrons and ions can be observed, for example, in the outer part of galaxy clusters 
(e.g., Markevitch et al. 1996; Fox & Loeb 1997; Ettori & Fabian 1998; Takizawa 1998). 

Dispersion relation (46) has a general form which permits us to investigate analytically 
different limiting cases. When the thermal pressure is larger than the magnetic pressure, 
the case which is satisfied in the intracluster medium, one can omit the first term on the 
right hand-side of Equation (46). Then, we can write this equation in the form 

uP_ _ R 

fcpf " (Teo + T,o)\/' ^ ' 

where Cg is the sound speed. It is followed from Equation (49) that in the case cu^ ^ kyC^ 
or R^ (Tgo + Tio) V, dispersion relation is given by 

y = 0. (50) 

In the opposite case, cu"^ <C kyC^ or <C (Tgo + T^q) V, we have 

R = Q. (51) 



9.1. SOLUTION OF DISPERSION RELATION V ^0 

When r2je,e« = 0, i.e. at the absence of energy exchange. Equation (50) gives 

{-iuj + + VLtc) {-ioo + ^Ti) = 0. (52) 

Thus, an instabihty can be generated by the electrons or ions. We see that Equation (52) 
has solutions which correspond to the isochoric ones in the MHD (Parker 1953; Field 1965). 
When fije,ei — > oo, the case of a strong energy coupling, we obtain 

_ • (^X + ^Te) ^ie + ^Ti {^ei + ^Tei) 
^ — ^ To o o \ ■ V / 
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Taking into account that = ^ei at Ueo — riio and ^^Tei — —^^ei {Teo — Tio) /2Teo 
(Braginskii 1965), this equation becomes the following: 

.2fl^ + 2flTe + ^T^{mo/TeO-l) 
1 + 3Tjo/TeO 

We see that if STio/T^o <S 1, then u — —i (20^ + 2QTe — ^Ti)- Thus, the ions can contribute 
to instability when Qtz > 0. In this case, condition of instability takes the form 



-eO ^^-i-iQ ll-o 

where we have only taken into account cooling functions of electrons and ions (see Section 2). 
We see from this condition that instability can also be possible in the electron temperature 
domain where dAe (Teo) /dT^o > 0. 



9.2. SOLUTION OF DISPERSION RELATION R^O 

We now consider the dispersion relation (51) which is appropriate in the case 
u;^ <^ kycl- This equation coincides with the dispersion relation in the fast sound speed 
regime (Nekrasov 2011). For reading convenience, we repeat here results given in the 
last paper. By using the temperature dependence of ~ Te^^'^, Equation (51) can be 
rewritten in the form, 

TeO {-ijOJ + 0.^ + Q.Te " ^ne) {-ioj + Qti + ^^ie) (55) 

+ Tio {-i'yu + Qti - ^ni) + + ^Te + 21^ie) 

3 

~ 2^*6 (^eO ~ Tio) 

The different limiting cases of Equation (55) arc given in subsections of 9.2. 



[-i'^U + ^Ti) 1 + 



iO 



eO 



0. 
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9.2.1. The case Qie^O 

If we do not take into account the energy exchange, = 0, and set Tgo = Tio, then we 
obtain equation 

270;^ + i [(7 + 1) (^^x + ^Te + ^Ti) - ^ne " ^ni] (56) 
— 2 {fl^ + Qtc) ^Ti + ^ne^Ti + ^ni (^x + ^Te) = 0. 

Neglecting the contribution of the ion coohng and heating, ^2^^ = flni = 0, we have 

^ = -77- [(7 + 1) (^X + ^Te) - ^ne] ■ 

It is easy to see that this solution is a mixture of isochoric and isobaric solutions (Parker 
1953; Field 1965) because we have taken into account the ion temperature perturbation. If 
we neglect the latter, i.e. neglect the second term ~ T^o in Equation (55), we obtain the 
usual isobaric solution 

UJ — —— {fl^ + flxe — ^ne) • 

We also see from Equation (55) that for short-wavelength perturbations when » 
cu, flxe, ^ne the thermal instability can arise due to the ion cooling function 

^ — "Tfi ; 7^ [(^eO + ^io) ^Ti — ^^iO^m] • 

9.2.2. The case flie oo 

When the frequency flie is much larger than other frequencies, 2Qie ^ ^, ^x' ^Te,i, and 
Teo — TiQ, then the dispersion relation becomes the following: 

a; = - — {fl^ + flTe - fine + f^Ti - f^ni) ■ (57) 



This is isobaric solution with the electron and ion coohng-heating. 

For different temperatures of electrons and ions, T^q ^ Tjo, we obtain 



• 7 
I— 
2 



Teo+(4 + 3^)T,o 



(3TeO — Tjo) {^x + ^Te) 



1 / ^0 

~ o i'^'^iO + ^eo) I ^ne + 7^;— ^ni 



In the case Teo 3> Tjo, this equation takes the form, 



In the opposite case, Teo ^ Tjo, we obtain the ion isobaric solution 



U! — {^Ti ~ ^m) 

7 



5.^.5. General case 



In a general case. Equation (55) can be written in the form 



goUJ + igiUJ - g2 = 0, 
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where 



^0 = 7(^60 + ^,0), 



r„„ + Tio I 4 + 3^ 



eO 



+ (STeO — Tjo) f^ie (^x + ^Te) ~ 



(60) 



1 / 

~ (2eO + STjo) ^ie I ^ne + ^niT^ 
^ \ J-eO 

Equation (59) can be solved numerically for known cooling-heating functions and 
temperatures. 



10. DISCUSSION 

Applying the operator V- to Equation (8), we can see that Ugi — rin, if the electron 
and ion number densities are only perturbed. However, in our general calculations in the 
Appendix, the values V ■ Vgi and V ■ Vji are considered as different in the cases Ueo 7^ tiio 
and — riio. Solving equations of motion, we do not need to use the condition nn — Ugi. 
These equations permit us to treat a general case nn ^ n^i that can be appropriate 
for multicomponent plasmas. To derive the dispersion relation, we find expressions for 
perturbed velocities and calculate the perturbed current which then is used in Equation 
(8) in his given form. The perturbed current does not contain density perturbations at the 
absent of background flow velocities. Equations (24)- (29) are justified under conditions 
(18) and (19) and can be appUcd for both riio = Ueo, rin = Uei and riio 7^ ^eo, nn ^ Uei 
(for multicomponent plasmas) cases. Taking into account condition (18), we obtain from 
Equations (24), (25), (28), and (29) that V • Va = V • Vgi = -{c/Bo){dEi^/dy) for 
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perturbations (44). 

From the results obtained above, we can estimate relative perturbations of number 
density and pressure in the fast dynamical or long- wavelength regime (18). Using, for 
example, Equations (2), (24) and (28) and keeping main terms, we find equation for the ion 
density perturbation, 

Hiodt Uci dy 

From Equation (21), it follows that 

P^l = {-K + l^^)—^^, (62) 

Uci dydt 

where we have used the relation mjFji = — rngFei. The value Pji is connected with the 
pressure perturbation pji as follows 

Pn = -— (63) 

mj-rijo at 

From Equations (61)-(63), we obtain 

driji ^ TjQ Pa 
riiodt {Xi - Hi) rrii pio ' 

Analogously, we have for electrons 

driel _ TeO Pel 

Uiodt {Xe- llejrriePeo' 

where rZei = Un. Thus, 

Pil Pel 



(64) 



(65) 



(Xi - Hi) rui (Ae - He) rUe ' 
Taking into account notations (23), we see that ni^eil^io ~ Pi,ei/Peo and pn/pio ~ Pei/Peo- 
Using the dispersion relation (44) in the case of neglect the magnetic field, we obtain from 
Equations (64) and (65) equation connecting the sum of pressures and density perturbation, 

d'rii, ^ d\ 
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It is easy to see from Equations (33), (38), and (44) that Equation (51) corresponds to 
isobaric regime where the sum of electron and ion pressure perturbations is smaller then 
electron or ion pressure perturbation (see also Nekrasov 2011). The perturbation of number 
density is due to electric drift (see Equation (61)). 

Dispersion relation (50) is satisfied in the case cu^ ^ ^y^l- Taking into account 
condition (18), unstable perturbations have ky < {mi/me) k^. Thus, the transverse 
wavelength A_l > {me/rriiY^'^ Xz and can be both less and larger than the longitudinal 
wavelength A^. From other side, in the case cu^ ^ ^ycf; dispersion relation (51) describes 
unstable perturbations strongly elongated along the magnetic field, ky ^ {mi/me) k^ or 
Xz ^ {mi/nie)^^^ X±. In this case, very thin filaments are generated. Thus, a wide spectrum 
of wavelengths of perturbations along and across the magnetic field can be formed in the 
framework of conditions (18), (19), (34), and (41). 

A general form of dispersion relation (46) including the thermal exchange and different 
temperatures allows us to consider various cases, which can be realized in real situations. 
We can investigate a weakly, strongly, and intermediate thermal coupling (see Sections 5.5 
and 5.6). In particular. Equations (52) and (56) are available for a weak thermal coupling, 
while Equations (53), (54), (57), and (58) are appropriate in the case of strong coupling. 
The intermediate case is described by Equation (59), where coefficients (60) contain both 
different temperatures and different cooling functions. 

We have shown that unstable perturbations have an electromagnetic nature (see 
Equation (44)). Thus, a consideration of only potential perturbations is in general not 
adequate. 



11. ASTROPHYSICAL IMPLICATIONS 
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We shortly outline some important points of our investigation for possible observations. 
We have found growth rates, which contain in the clear form the separate contribution 
of cooling functions of electrons and ions. It is obvious that both components (in 
multicomponent media also dust grains, neutrals, and so on) can result in thermal 
instability in the same extent. This fact considerably extends possibilities for the medium 
to become unstable. In this connection, it is important to know the functional dependence 
of cooling functions on the temperature and density for each species. Unfortunately, at 
present, there is not sufficient information on this subject in astrophysical literature. 
The range of scale lengths of unstable perturbations can enlarge due to contribution to 
instability of other species except electrons. For example, short-wavelength perturbations, 
which must be stable because of a large electron thermal conduction, can be unstable 
due to contribution of ions to coohng of medium (see Section 9). In the long- wavelength 
regime (18), scale sizes of unstable perturbations across the magnetic field can have a wide 
spectrum and be, in particular, very elongated along the magnetic field. Such filaments 
are observed in galaxy clusters (e.g., Conselice et al. 2001; Salome et al. 2006) and in the 
solar corona (e.g., Tandberg-Hanssen 1974; Karpen et al. 1989). Different temperatures of 
electrons and ions assumed in this paper can be observed in galaxy clusters (Markevitch 
et al. 1996; Fox & Loeb 1997; Ettori & Fabian 1998; Takizawa 1998). This fact proves 
that dynamical and statistical processes could have timescales of the same order. It is clear 
that real situations in astrophysical objects are much more complicated to be captured 
by simplified theories. Knowledge of fundamental processes and more detailed conditions 
from observations are very important for theoretical models and, in particular, for further 
investigation of thermal instabilities. 



12. CONCLUSION 
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We have studied thermal instabihty in the electron-ion magnetized plasma which is 
relevant to galaxy clusters, solar corona, and other two-component astrophysical objects. 
The multicomponent plasma approach have been applied to derive the dispersion relation 
for the condensation mode in the case in which the dynamical time is smaller than a time 
the particles need to cover the wavelength of perturbations along the magnetic field due 
to their thermal velocity. Our dispersion relation takes into account the electron and ion 
cooling-heating functions, collisions in momentum equations, energy exchange in thermal 
equations, different background temperatures of electrons and ions, and perturbation of 
energy exchange collision frequency due to density and temperature perturbations. Different 
limiting cases of dispersion relation have been considered and simple expressions for growth 
rates have been obtained. We have shown that perturbations have an electromagnetic 
nature. We have found that at conditions under consideration transverse scale sizes of 
unstable perturbations can have a wide spectrum relatively to longitudinal scale sizes 
and, in particular, form very thin filaments. General expressions for dynamical variables 
obtained in this paper can be applied for astrophysical and laboratory plasmas also 
containing the neutrals, dust grains, and other species. The results obtained can be useful 
for interpretation of observations of dense cold regions in astrophysical objects. 

In this paper, we have investigated the linear stage of thermal instability in the 
multicomponent medium. The instability development can result in plasma turbulence 
when transport coefficients become dependent not on Coulomb collisions but on the energy 
of turbulence. In this nonlinear consideration of a problem is necessary. 
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A. APPENDIX 
A.l. Perturbed velocities of species 

In the linear approximation, Equation (1) for the perturbed velocity Vji takes the form 

^ = _^+F,,+^v,,xB„, (Al) 

ot mjUjo rrijC 

where pji — rijoTji + rijiTjo. Prom this equation, we can find solutions for the components 
of Vji. For simplicity, we assume that d/dx = because a system is symmetric in the 
transverse direction relative to the ^-axis. Then, the x-component of Equation (Al) gives 



Fjlx+l^cjVjly, (A2) 



dt 

where u^-j = qjEo/mjC is the cyclotron frequency. Differentiating Equation (Al) over t and 
using Equation (2) in the linear approximation and Equations (15), (16), and (A2), we 
obtain for the y-component of Equation (Al) 



where 



^ + ^cjj vjiy = + Qjiy^ (A3) 



p G2 ^ j '^eo Gi 1 ^ f A4) 

Drrip dt V mg Drriedt ' 



Drrii dt \ rrii Drrii dt 
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The value Pji is connected with the pressure perturbation (see Equation (Al)). Using 
Equations (A2) and (A3), we find 



d 



Ucjdt 



dy 



(A5) 



In Equations (A3) and (A5), notations 



Qjly — ^cjFjlx ~l~ 



dF, 



dFjix 



(A6) 



dt 



are introduced. We see from these equations that the thermal pressure efi^ect on the velocity 
Viix is much larger than that on v^y when d/dt <^ uici- The z-component of Equation (Al) 
can be written in the form 



dz 



+ 



dt 



(AT) 



A. 2. Calculation of V • Vji and Pji 

We have 



V-v,^r^ + ^. (A8) 



Using Equations (A3), (A4), (A7), and (A8), we obtain 



Here 



Lie V • Vei + L2eV • Va = H^l, (A9) 

LiiV ■ Vil + L2iV ■ Vei = Ha- 
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and operators Lij and L2j are the following: 



Lie 



92 



Gi d 
Drrie dt 
Ga d 



TeO 

rrii Drrii dt 



^3e 



G2 d 

Dnie dt 



1 J^2i 



+ 



^3 d 

^* Drrii dt ' 
92 



Prom a system of equations ( A9) , we find 



^ , 2 

^2 +^cj / Q^2 



LV • Vel = —L^eHii + LiiHei, 

LV ■ Vji = —L2iHei + Lie -f^a, 



(All) 



(A12) 



where 

L = L\eL\i — L2eL2i. (-^13) 

The values Pei and Pn can be found, substituting solutions (A12) into expressions (A4), 



LPa 



Drrii dt \ 
Drrii dt * \ m,- 



Gj d 
Drrii dt 
_G4_d_ 
Drrii dt 



L2i 



Hii 

Hel- 



(A14) 



d ^ /TeO 
Drrie dt ^ V 
G2 d fTeo 



Drrie dt 



Gi d 
Drrie dt 

Gi d 
Drrie dt 



Lii 



■i2e 



Hil, 



(A15) 



A. 3. Expressions for D and Gi,2,3,4 

We now give expressions for values defined by a system (17): 

/ d \ d \ d'^ ( d \ cP' 

D=[-^ + ^X + ^Tej [-Ql + ^Ti + nie) + {^ei + ^Tei) + ^Ti\ ^, (A16) 
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Gi = Teo Vine - (7 - 1 



d_ 
di 



^^+^T^ + ^ + (Teo - T,o) + ^Ti) ^, (A17) 



(-4 



di 



- + Qei {Teo - Tio) f ^ + flT^] ^, (AI8) 
— - (Teo - T,o) ( — + + l^Te j ^, (Al9) 



( d 

— Qje (2^eO — 2^io) f ^ + + ^Te j 



^^m - (7 - 1) 



d_ 
di 



d_ 
di 



(A20) 



A.4. Simplification of Equations (A14) and (A15) 



We further calculate coefficients by Hji in Equations (A14) and (A15). Using 
expressions (All), we find 



Gs d ^ , fT,o 



Drrii dt 



Ga d \ ^ G^ f d^ o\ 93 



m,- Drrii dt 



Dm, \dt^ '^'^"'1 dt^ 



(A21) 



and 

Gs d 
Drrii dt 



L2e + 



iO 



Ga d 



rrii Drrii dt 



Lie 



D 



Tio GAd\ ( d 



m,- m,- dt J \dt 



d^ 



+ 



dt"^ DmeTrii 



LseK. 

(A22) 



In Equation (A22), we have introduced notation 

1 d 

-f^ = {G2G3 - GiGa) — + (reoG4 + TioGi) — - DTeoTiQ. 



(A23) 



Calculations show that the value {G2G3 — G1G4) has a simple form, i.e.. 



— {G2G2, — G1G4) — ^ie {TgQ — Tjo) TeO 



d 

^ne - (7 - 1) ^ 



+ ^ei {TjQ — Teo) Tj 



d 

(A24) 



~ TeoTio 



^ne - (7 - 1) 



d_' 
di 



^ni - (7 - 1) 



d_' 
di 
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Using expressions (A16), (A17), (A20), and (A24), we obtain for the operator K (A23) the 
simple form, 



92 



92 



where notations 



d 
d 



(A25) 



(A26) 



are introduced. Using Equations (A21) and (A22), Equation (A14) for Pa takes the form, 



DLPa 



D- 



iO 



Ga d 



rrii 



G. ( & 



rriidt J \dt^ 



+ 



LseK 



Hii 



Analogous consideration of Equation (A15) leads to the following equation for Pei. 



DLPei 



rUe rriedt ) \ dt"^ dt'^ rriime 



Hel 



G, ra' 



^2 I u ,^o\d^ 



Operators 



rrie me dt ' 
rrij rrii dt 



can be found by using Equations (A16), (A17), (A20), and (A26) 



■_eO 



G^d_^T^, 
nif, dt nip 



rrif. \ot I at^ 



- = -^W, ( ^ + + Ore + ^^ei + ^Tei 



.Tio G4 d _ TiQ 
m 



d_ 

at 

d_ 
dt 



dt^ 



+ ^Qie ( 7^7 + + ^Te 1 — • 



dt^' 



(A27) 



(A28) 



(A29) 
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A.5. Operator L in a general form 

Using expressions (All), we find from Equation (A13) 

L^M-N ^ L^^L^,K, (A30) 

mefniD 

where 

V me DrUe dtJ^\dt^^'^''J dt^ V mi Drm dt J ' 
and K is defined by Equation (A25). 
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